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The polarization state of the optical electromagnetic field lying several nanometers above complex 
dielectric-air interfaces reveals the intricate light -matter interaction that occurs in the near-field 
zone. From the experimental point of view, access to this information is not direct and can only 
be extracted from an analysis of the polarization state of the detected light. These polarization 
states can be calculated by different numerical methods well-suited to near-field optics. In this 
paper, we apply two different techniques (Localized Green Function Method and Differential Theory 
of Gratings) to separate each polarisation component associated with both electric and magnetic 
optical near-fields produced by nanometer sized objects. A simple dipolar model is used to achieve 
insight into the physical origin of the near-field polarization state. In a second stage, accurate 
numerical simulations of field maps complete data produced by analytical models. We conclude this 
study by demonstrating the role played by the near-field polarization in the formation of the local 
density of states. 

PACS numbers: 42.79.G, 42.82.E, 07.79.F 



I. INTRODUCTION 



Light interactions with dielectric or metallic surfaces 
displaying well-defined subwavelength-sized structures 
(natural or lithographically designed) give rise to unusual 

optical effects@, I> & & 6 I' I E E & Thc 

recently observed "light confinement state" in which the 
light field is trapped by individual surface defects, be- 
longs to this class of phenomena^. Although, with usual 
dielectric materials (silica for example), the local near- 
field intensity variations observed around the particles (or 
structures) remains moderate over the optical spectrum 
(between 10 to 40 per cent of the incident light inten- 
sity), these variations can nevertheless be easily mapped 
with the tip of a Photon Scanning Tunneling Microscope 
(PSTM)Q. The images recorded with this technique re- 
veal dramatic changes when passing from the TM (trans- 
verse magnetic) to the TE (transverse electric) -polarized 
modes. In general, TM-polarized light tends to display 
larger contrast than TE-polarized light. ^From the ex- 
perimental point of view, the definition of thc polar- 
ization direction of the incident and detected intensity 
must be defined with respect to a unique incident plane. 
About six years ago Van Hulst and collaborators pro- 
posed a clever probe configuration devoted to polariza- 
tion mapping [ fi3| . These authors performed these mea- 
surements by using a combined PSTM/AFM microscope 
in which detection is implemented by a microfabricated 
silicon-nitride probe. ^From this technique, polarization 
contrast is extracted by changing the polarization direc- 
tions of both the incident and the detected light. The 



main findings gathered in this work]13|], concern the rel- 
ative efficiency of the four excitation-detection possibil- 
ities (TE/TE,TE/TM,TM/TE,and TM/TM) to record 
a highly resolved PSTM image. In particular, the ef- 
ficiency of the TM/TM acquisition mode is well repro- 
duced. Although a complete interpretation of this work 
requires a realistic numerical implementation of the com- 
bined AFM/PSTM probe tip, we can obtain useful in- 
formation by analyzing the near-field polarization state 
versus the polarization state of the illumination mode. 

In addition, in closely related contexts, the control 
of the near-field polarization state provides an interest- 
ing and versatile tool for generating powerful applica- 
tions (tunneling time measurements |14 ] , highly-resolved 
microscopy and spectroscopy Q, surface plasmon reso- 
nance spectroscopy of molecular adlayers 
11 " 



tics 



20 



If 



atom op- 
More precisely, in the field of atom optics 
and intcrfcrometry, one is interested in building diffrac- 
tion gratings that can play the role of beam splitters. 
Several devices have been successfully realised, ranging 
from mechanical transmission gratings to light standing 
waves in free space or evanescent for a prism. For a 
general review the reader is referred to |16| . To circum- 
vent some theoritical limitations |l7| it has been recently 
proposed [fbsf to use micrometer sized metallic stripes 
to shape the evanescent field. A full near-field, metal- 
lic/dielectric approach open obviously new perspectives. 
In particular the spacing period is no longer linked to the 
atomic optical transition and a reflection stucture cannot 
clug. More, higher harmonics in the optical evanescent 
field can be tuned to produce for example a blazed atomic 
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grating. Nevertheless, the optical potential is strongly 
related to the light field polarisation, which is a farther 
motivation for the present study. We will begin our the- 
oretical analysis with a simple dipolar scheme in which 
the main experimental parameters (incident angle, opti- 
cal index, polarization of the incident light, ...) appear 
explicitly pl| |22] , |23|. In a second stage, these results will 
be completed with an ab-initio approach allowing objects 
of arbitrary shape to be treated exactly. 



II. POLARIZATION OF THE LIGHT ABOVE A 
SINGLE DIELECTRIC PARTICLE 

To illustrate the coupling between a polarized incident 
wave and a small spherical object lying on the sample, 
we consider the simple dipolar model depicted in Fig. 1. 
The substrate modifies the polarizability ao(u>) of the 
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FIG. 1: Schematic drawing of the model system used in the 
present section. A transparent substrate, of optical index n 
= ei supports a small dielectric sphere of radius R. The sys- 
tem is illuminated in total internal reflection with an incident 
angle 9q and R = (0,0,7?). xOz is the incident plane. 



particle. We have then 

a eff {R,uj) =a {uj) •M(R,w) 

with 

M(R, oj) = [1 - S S (R, R, u) • aoMr 1 



(1) 



(2) 



where S S (R, R, u>) is the nonretarded propagator associ- 
ated with the bare surface, and R = (0, 0, R) labels the 
particle location. Within this description, the optical 
properties of the spherical particle-surface supersystem 
is described in terms of a "dressed" polarizability |2^, 
The analytical form of a e ^ can be derived from Eq. (10) 
of reference 5. This dyadic tensor remains diagonal with 
two independent components a e /^ (perpendicular to the 
interface) and af/* (parallel to the interface): 



eff 



8{n 2 + l)a (uj)R 3 



8(n 2 + 1)R 3 -a (uj)(n 2 - 1) 



(3) 



(R,w) 



4(n 2 + l)a (io)R 3 



(4) 



4(n 2 + l) J R 3 -a (w)(n 2 -l) 
where n is the optical index of refraction of the substrate. 

A. New field components in the near— field domain 

At an observation point r located above the sample 
[z > 0) and in the immediate proximity of the particle, 
the incident light field is locally distorted. As illustrated 
in Ref. p5| , these distortions generate not only a pro- 
found modification of the intensity level (both electric 
and magnetic), but also a complete change of the po- 
larization state. At subwavelength distances from the 
scatterers, we expect therefore to observe the occurrence 
of new components that were absent in the incident field 
{Eo(r, t); Bo(r, t)}. The physical origin of this polariza- 
tion transfer can be easily understood if we introduce 
the two relevant field propagators So and Qo that es- 
tablish the physical link between the oscillating dipole 
fi(t) = a e ^(R, wo) • E (R, t) and the new near-field state 
{E(r, t); B(r, t)} generated above the particle 



E(r, t) = E (r, f) + S (r, R) • a e ^(R, uj ) ■ E (R, t) 



(5) 



and 



B(r, t) = B (r, t) + Q„(r, R, w ) ■ a eff (R, lu ) ■ E (R, t) 

(6) 

where in the near-field zone, i.e. when |r — R| < A = 
2irc/uj (h 

S„( r .R) ^' r - R » f - R »-"- R ' 21 (7) 



|r-R| 5 



and 

Q (r,R, w ) 
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(8) 



The discussion of the two equations (5) and (6) can be 
made easier when considering specific examples. 



B. Illumination with a TE— polarized surface wave 

In this incident polarization mode, the electric field is 
directed along the (OY) axis. We have then Eo(R, t) = 
(0, Eo(t), 0). Let us see what happens when the observa- 
tion point moves along the diagonal straigth line (A-B) 
schematized in figure 2. For a given (3 angle, the intro- 
duction of the position vector r = (r cos(/3), r sin(/3), z) 
(where r varies between — oo and oo) along the line (A- 
B), allows the electric field polarization evolution to be 
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FIG. 2: Top view of the single spherical system depicted in 
figure jjj. The straight line (A-B) that passes over the sphere 
center at a constant height z is used to evaluate the polariza- 
tion change produced by the particle. 



observed when passing over the particle. Each cartesian 
component can be simply extracted from Eq. (5). This 
leads to three analytical relations: 



E x it) = E (t) 



af/ f (R,w)3r 2 cos(/3) sin(/3) 



\r 2 + (z - R) 2 \ 5 / 2 
E y (t)=E (t){l + al ff (R,w )T(f3)}, 



with 



n/3) 



and 



E z (t) = Eo(t)- 



2r 2 sin 2 {P) - r 2 cos 2 (P) - (z - R) 2 
\r 2 + (z - R) 2 \ 5 / 2 

iu // (R,tJo)3rsin(/3)(2 - R) 



\r 2 + (z - R) 2 \ 5 / 2 



(9) 
(10) 

(11) 
(12) 



Some interesting features can be deduced from these 
equations, (i) First, we observe that the polarisation is 
not modified when the observation point is perpendicular 
the particle center (i.e. when r = 0). The set of equations 
reduces then to 



E x (t) = 0, 



Ey(t) = E t 



oWi 1 r Biz / 



z - R[ 



and 



E z (t)=0. 



(13) 



(14) 



(15) 



Clearly the effective polarizability reduces the electric 
field magnitude compared to its initial value. This fact 
explains the observation of contrasted dark zones above 
small particles lighted with the TE polarization (c.f. Fig. 
(4)). (ii) Around the particle (when R/2 < r < 2R) two 
new components, namely E x (t) and E z (t), define a new 
local polarization state. These components vanish again 
when the observation point moves away from the particle. 
Similar relations can be derived for the TM-polarized 
mode from Eq. (5). 

To conclude this section, let us examine what happens 
with the magnetic field part (cf. Eq. 6). Since in the 
reference system of Fig. (1), the incident magnetic field 
displays two components different from zero, Bo(r,t) = 
{Box(t), 0, Bo* CO) j w e can write 



B x (t) — B 0x (t) 



iu (z - R) eff 



c[r 2 + {z- R) 2 fl 2 



and 



B z (t) = B 0z (t) 



By{t) = 0, 



itoor cos(P) 



al JJ E (t), (16) 



(17) 



i? )2]3/2 Q! | 



eff 



E (t). (18) 



Unlike what happens with the electric field, the particle 
does not produce new magnetic field components in the 
near-field. In this case, the polarization change corre- 
sponds to a different balance in the initial components. 
It is important to recall that we have used the dipole ap- 
proximation to describe the particle-field interaction, i.e. 
the size of the particle is assumed to be small compared 
to the wavelength of light. In a more realistic calcula- 
tion, with nanostructures of characteristic dimension ss 
100 nanometers this result is not rigorously exact. How- 
ever, we still expect, in TE polarization mode, a negli- 
gible particle contribution to B y compared to the total 
magnetic field intensity. 



III. AB INITIO STUDY OF THE NEAR-FIELD 
POLARIZATION STATE 

Analytical results presented in the previous section 
supply qualitative information about the spatial polar- 
ization state distribution. In a recent paper, analysis of 
polarization effects was proposed in the context of near- 
field optics in which a limited number of single parti- 
cles were investigated |2r| . Since in many practical situ- 
ations experimentalists are interested in lithographically 
designed structures, these preliminary analyses must be 
completed by ab initio procedures for solving Maxwell's 
equations. 



A. Localized objects 

Recently, theoretical modelling in the vicinity of lo- 
calized objects was performed in the framework of the 
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field susceptibility method|27], |2f|. Today, this method 
is one of the most versatile and reliable numerical tech- 
niques to solve the full set of Maxwell equations for 
the typical parameters of near-field optics. It works 
well even for metallic nanostructures (see for example 
references |t], ||, [To)). This approach (called the Direct 
Space Integral Equation Method (DSIEM)) is based on 
the knowledge of the retarded dyadic tensor S(r, r', oS) as- 
sociated with a reference svstem which, in our problem, 

. The numerical procedure 



is a flat silica surface) 29 



considers any object deposited on the surface as a lo- 
calized perturbation which is discretized in direct space 
over a predefined volume mesh of N points {Ri}. In a 
first step, the electric field distribution E(Ri,w) is de- 
termined self-consistently inside the perturbations (i.e., 
the source field). At this stage, a renormalization proce- 
dure associated to the depolarization effect is applied to 
take care of the self-interaction of each discretization cell. 
The final step relies on the Huygens-Fresnel principle to 
compute the electromagnetic field E(r, uj) on the basis of 
the knowledge of the field inside the localized perturba- 
tions E(Ri,w). The two main computational steps can 
be summarized as follows: 



In the numerical work to be discussed in this section, 
the retarded propagators S and Q have been chosen in 



reference! 31 



The test-object we consider in this section is the word 
OPTICS engraved at the surface of a TM2 layer de- 
posited on a silica surface. Intentionally we have chosen 
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(i) Local field computation inside the source field 

E(R i ,w) = 5^/C(R i ,Rj,w)-E (R J ,w), (19) 

3 

where /C labels the generalized field propagator of the 
entire system (localized object plus bare silica surface). 
In the {R^ Rj} representation it is given by 

^(Ri.Rj.w) = dij +v j S(Ri,TL j ,u) ■ x(R^), (20) 

where \ represents the electric susceptibility of the lo- 
calized object, Vj is the volume of the elementary dis- 
cretization cell, and S is the field-susceptibility of the 
entire system. This last quantity is usually computed by 
solving Dyson's equation: 

5(Ri,R J; w) = S(R i ,R j ,w) + (21) 

^i>fcS(Ri,R fe ,w) • x(Rfc,w) • S(Rfc,R 3 ,w), 

fe 

(ii) Electric and magnetic near-field mapping computa- 
tion around the source field region 

E(r,w) =E (r,w)+ (22) 
«<S(r, Ri, w) • x(Ri, w) • E(Rj, w). 

i 

and 

B(r,w) =B (r,cj)+ (23) 
^?; 4 Q(r,R 4 ,w) • x(Ri,w) • E(Rj,w). 

i 



200 400 

FIG. 3: Top view of the pattern used in the simulations of 
the ab initio studies. The computational window is 500 x 500 
nm 2 . 

a planar structure devoid of any symmetry. In order to 
gain more insight in the polarization changes occurring 
around such complex lithographically designed nanos- 
tructures, we analyze in Figs. (4)- (7) the electric and 
magnetic near-field intensities generated by each carte- 
sian component (E 2 X , E 2 ) and (B%, B 2 , B 2 Z ). For 
comparison, the square modulii are also provided. 



1. Dielectric materials 

The lateral dimensions of the object are given in Fig. 
(^) . The thickness and the optical index of the TiOi pat- 
tern is 20 nm and 2.1, respectively. The wavelength of 
the incident laser is 633 nm. All fields are computed 10 
nm above the surface of the structure, i.e. 30 nm above 
the glass-air interface. The incident light is a TM/TE 
polarized evanescent surface wave traveling along the Ox 
axis. This illumination condition is used in the Photon 
Scanning Tunneling Microscope (PSTM). Some general 
comments can be made about these results. First, all 
components of both the electric and magnetic fields have 
been excited in the near-zone. The occurence of these 
new components is a pure near-field effect because it is 
always localized around the structures. In Figs. (|J) and 
we display the electric field part. As predicted in sec- 
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polarized surface wave. In agreement with the PSTM re- 
sults, numerous regions appear with a dark contrast and 
a moderate intensity level. 

As expected, the excellent image-object relation cur- 
rently observed in the TM-polarized mode is mainly pro- 
vided by the field component E z normal to the object. 
The two other contributions tend to slightly degrade the 
total pattern E 2 composed by the superposition of the 
three maps E 2 ., E 2 and E 2 . 

The magnetic near-field intensity maps (cf. Figs. 6, 7) 
also show a significant confinement of the magnetic field 
over the particle which reverses the contrast with respect 
to the electric map. Similarly to what happens with the 
electric field, the role played by the additional compo- 
nents can degrade the topographic information contained 
in the complete field maps. Notice in Fig. (6) that, as 
mentioned in section II, the new y-component of the mag- 
netic field is very small compared to the total magnetic 
field. 



FIG. 4: Grey scale representation of the electric field dis- 
tribution computed above the topographical object depicted 
in Fig. 3. The calculation is performed in the TE-polarized 
mode and the arrow indicates the propagation direction of the 
surface wave. Extreme values of the components of the elec- 
tric field (normalized by the incident field) are : 0.000 (min) 
and 0.153E-1 (max) for El, 0.674 and 1.67 for E%, 0.000 and 
0.747E-1 for El, 0.681 and 1.68 for E 2 . 
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FIG. 6: Grey scale representation of the magnetic field dis- 
tribution computed above the topographical object depicted 
in figure 3. The calculation is performed in the TE-polarized 
mode. Each map is normalized with respect to the incident 
magnetic field intensity. Extreme values are : 0.287 (min) and 
0.396 (max) for Bl, 0.696E-13 and 0.490E-4 for B 2 , 0.684 and 
0.772 for B 2 Z , 0.973 and 1.145 for B 2 . 



FIG. 5: Same as figure 4, but in TM-polarized mode. Ex- 
treme values are : 0.206 (min) and 0.436 (max) for El, 
0.700E-9 and 0.216E-1 for E 2 , 0.613 and 1.06 for E 2 Z , 0.908 
and 1.31 for E 2 . 



tion 2, we recover the appearance of two additional com- 
ponents, E x and E z , when the object is excited by a TE- 



2. Metallic materials 

In the above formalism, the only parameter distin- 
guishing metallic from dielectric objects is the linear 
susceptibility \{y',uo). Alternative procedures can be 
adopted to describe the metallic susceptibility. For ex- 
ample, a direct route would consist in expanding the sus- 
ceptibility in a multipolar series around the geometrical 
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FIG. 7: Same as figure 6, but in TM-polarized mode. Ex- 
treme values are : 0.109E-7 (min) and 0.856E-2 (max) for B 2 , 
0.710 and 1.11 for B 2 , 0.107E-7 and 0.190E-1 for B 2 , 0.711 
and 1.11 for B 2 . 

center of the metallic particle. While this scheme allows 
non-local and quantum size effects to be included, it is 
nevertheless restricted to simple particle shapes (spheres, 
ellOipsoids, etc.). When dealing with spherical metallic 
clusters having a typical radius below 15 nm this de- 
scription is mandatory and can be easily included in the 
DSIEM formalism |32|. For the applications discussed in 
this paper, involving lithographically designed metallic 
structures larger than this critical size, we can adopt the 
discretization of x( r 'i w ) over all the volume occupied by 
the particle. In this case, the local susceptibility is just 
related to the metal optical index n by the relation Q: 



In the visible range, the numerical data for describing 
both real and imaginary parts of n have been tabulated 
by Palik||| for different metals. We present in Fig. (||) 
a gray scale representation of the electric field distribu- 
tion computed above the topographical object depicted 
in Fig. 3. In this case, the high optical metal index gen- 
erates complex field patterns without clear relation to 
the topography. Furthermore the possible excitation of 
localized plasmons reinforces this phenomena and some 
parts of the localized metal pattern (e.g. the corners) can 
even behave as an efficient light sources. 

B. Periodic surface structures 

When working with periodic surface structures, the lo- 
calized Green's function method described above is no 



FIG. 8: Grey scale representation of the electric field distri- 
bution computed above the topographical object depicted in 
Fig. 3. In this application the supported structure is metallic 
(gold). The incident wavelength is 760 nm. The calculation 
is performed in the plane zq = 30 nm in the TM-polarized 
mode and the arrow indicates the propagation direction of 
the surface wave. Extremal values of the components of the 
normalized electric field are : 0.118E-2 (min) and 3.13 (max) 
for El, 0.146E-5 and 1.72 for E 2 , 0.237 and 10.1 for E 2 , 0.851 
and 10.6 for E 2 . 

longer applicable. But any figure can be decomposed in 
points (direct space) or Fourier components ( reciproqual 
space) . We thus can use two methods wether the real ob- 
ject contain few points or few Fourier components. In the 
case of periodic gratings, the field distribution can be in- 
vestigated with this second class of methods |m], ^5[ [36| . 
The so-called differential theory of gratings (DTG) was 
originally developed twenty years ago to predict the effi- 
ciencies of one- and two-dimensional diffraction gratings. 
Based on a rigorous treatment of Maxwell's equations, 
this method can also be used efficiently to determine the 
optical near-field scattered by three dimensional periodic 
objects. In the following subsection, in order to avoid a 
complete presentation of this well established technique, 
we will only summarize the essential steps of the compu- 
tational procedure. 

As in previous sections, we are interested in the elec- 
tromagnetic near-field diffracted above objects engraved 
on an interface illuminated by total internal reflection. 
When using the DTG method Q], the electromagnetic 
field above the grating can be expanded in a Fourier se- 
ries 

+OC +OG 

A W = J2 J2 Ap^e 11 "^^"" 1 , (25) 

p— — 00 q— — OG 

where r = (1, z) = (x,y,z), A(r) represents either the 
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electric field E(r) or the magnetic field B(r). The 3D- 
wave vectors k p q — (k|| pg , 7 p . g ) , associated with the har- 
monic (p, q) obey the well-known dispersion equation 



2; 2 



(26) 



The set of wave vector k|| p A parallel to the surface are 
simply defined for each couple of integer numbers (p, q) 
by 



k| = j nk 0x +P-^- ) n. 



2tt . 

nk 0y + q-y- ] u 



l y > 



(27) 



where d x and d y denote respectively the period of the 



grating along the Ox- and Oy-directions. From Eq. (26) 



it may be seen that the coefficient 7 p . g may be either 
real or purely imaginary. Real values of 7 P . g correspond 
to radiative harmonics while imaginary values introduce 
evanescent components in the expansion (pq). 

In a general way, the six components of the electromag- 
netic field A(r) can be deduced from two independent pa- 
rameters usually named the principal components. Let us 
choose, for example, the y-components E y (r) and B y (r) 
as principal components. It is a simple matter to show 
that the Fourier y-components of the field just above the 
surface of objects can be expressed as a linear combina- 
tion of the y-components of the incident field: 



yp,q 
yp,q 



= T 



EE 



Oy 
Oy 



n-EB T3 

n-BB R 
J pq 



Oy 



(28) 



The transmission coefficients T EE , T EB , T BE and T BB 
describe the coupling between the electric and magnetic 
harmonics composing the scattered and the incident field. 
These coefficients depend both on the geometry of the 
sample and on the angular conditions of incidence but 
not on the polarization of the incident light. The po- 
larization of the incident plane wave is controlled by the 
values of Bq v and Eo y ■ From a numerical point of view, 
the transmission coefficients are obtained by the inver- 
sion of a complex square matrix whose dimension is 2Nt 
x 2Nt (where Nt is the total number of harmonics used 
to describe the scattered field in Eq. (pq)). Columns 
of this matrix contain the Fourier y-components of the 
electromagnetic field which would have illuminated the 
periodical objects in order to obtain a pure harmonic 
field (p, q) just above the nanostructure. A detailed de- 
scription of the calculation of the matrix elements can be 
found in Refs. 
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|3J. With Eq. (|28|), we can calculate 
all the Fourier components of electric and magnetic fields 
just above the objects, which are used as initial condi- 
tions to obtain the field anywhere. General remarks 
about contrast, relative intensities and image-o bjec t re- 
lation have been made previously (see sections |ll B| and 
IIIA1), therefore we only highlight the specific electro- 
magnetic properties of periodic structures. 

We studied a lattice of 100 x 100 x 100 nm 3 Ti0 2 dots, 
separated by 150 nanometers. A strong localization of 
the electric field appears above the pads in TM-polarized 
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FIG. 9: Grey scale representation of the electric field dis- 
tribution computed above a regular lattice of square shaped 
dielectric pads. The calculation is performed in the plane zq 
= 50 nm in the TE-polarized mode, and the arrow indicates 
the propagation direction of the surface wave. Extremal val- 
ues of the components of the normalized electric field are : 
0.43E-8 (min) and 0.875E-2 (max) for E%, 0.266 and 0.713 
for El, 0.141E-7 and 0.644E-1 for El, 0.275 and 0.741 for E 2 . 



x dim) 



x dim) 




y dim) 



y dim) 



FIG. 10: Same as figure (9), but in the TM-polarized mode. 
Extremal values of the components of the normalized electric 
field are : 0.153 (min) and 0.454 (max) for E 2 X , 0.26E-8 and 
0.458E-1 for E 2 , 0.513 and 1.49 for E 2 , 0.802 and 1.58 for E 2 . 



mode or between the pads in TE-polarized mode. More- 
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over, a careful analysis of the different components shows 
that it is possible to create a particular field map such 
as field lines oriented along y-axis (E y in TE-mode) or 
x-axis (E x in TM-mode) , or periodic field spots with very 
different characteristics (spot size, periodicity and shape) 
considering the other components in the two polarization 
modes. These particular field components distribution 
could have a great interest for the interaction of cold 
atoms with optical evanescent waves p7L pq] . 



IV. LOCAL DENSITY OF STATE AND 
POLARIZATION EFFECTS 

Unlike what happens with electronic surface states, the 
Local Density of photonic States (the so-called photonic 
LDOS) contains information related to the polarization 
of the excitation field. It is well established that the 
density of states near surfaces plays a significant role in 
near-field physics In particular, the photonic LDOS 
is a useful concept for the inter pre tation of fluorescence 
decay rate in the very near-field |p9| and could help in un- 
derstanding image formation produced by illuminating- 
probe SNOM pQ| . By referring to the electric field, we 
deduce the optical LDOS from the field-susceptibility of 
the entire system (pl ane s urface plus supported nanos- 
tructures, see section [II A ) |4(i| , |4l|| 



p(r,u) 



1 



2tt 2 w 



Im [TrS(r, r, us) 



(29) 



In this expression, the optical LDOS is related to the 
square modulus of the electric field associated with all 
electromagnetic eigenmodes of angular frequency u>. Be- 
cause of the vectorial character of electromagnetic fields, 
it is very useful to introduce three polarized optical LDOS 
also called partial LDOS, so that (hi 



p(r,w) = p x (r,w) + p y (r,u) + p z (r,uj), (30) 
I 



Pi{r,Lu) = T ^-ImS ii (r,r,uj), i = x,y,z(31) 



The three different polarized optical LDOS computed 
over a pattern made of three dielectric cylinders of optical 
index 2.1 are represented on figure [H]. At this stage, we 
can made some general remark. In the configuration in- 
vestigated in figure ([n]) we can verify that both X-LDOS 
and Y-LDOS are almost identical with just a rotation of 
7r/2 between them. Nevertheless, this properties does 
not subsist any more if we reinforce the optical coupling 
between the dielectric posts. More precisely, certain ar- 
rangements of nanoscale pillars (circular or ellipsoidal) 
can force light waves into states generated by the col- 
lective coupling between the pillars leading to significant 
difference between X-LDOS and Y-LDOS maps. 



Moreover, although optical LDOS characterizes the 
spectroscopic properties of an electromagnetic system in- 
dependently of the illumination mode, it is interesting to 
note the strong analogy between the dark contrast over 




FIG. 11: a) Top view of the pattern used in the simulations 
of section |Vj consisting of three cylinders of 100 nm diameter 
and 60 nm height, b) z-ldos, c) x-ldos and d) y-ldos 150 nm 
above the surface at the wavelength A = 2n/u = 630 nm. 



the pads obtained for Y-LDOS and the dark contrast 
that appears in the case of s-polarized illumination mode 
(i-e inc ident electric field along the y-axis) observed in 
section 



[IB 



In addition, as in the case of optical near- 
field maps discussed in the previous section, the polarized 
LDOS's can display significant discrepancies relatively to 
the shapes of the underlying objects. 

Finally, let us note that partial LDOS's are not only 
a powerful mathematical tool but can easily be linked to 
the physical properties of electromagnetic systems. The 
most famous example is the fluorescence lifetime r of a 
molecule near an interface which critically depends on 
the spatial LDOS variation GJ, El EJ 



V. CONCLUSION 

On the basis of both simple analytical model and so- 
phisticated 3D Maxwell's equations solvers this paper has 
focussed on the unusual behaviour of the light polariza- 
tion in the near-field. When subwavelength patterned 
objects are excited by a surface wave of well-defined po- 
larization, a complex rotation of the light polarization 
state can be expected in the near zone. This phenomenon 
localized around the scatterer is a typical near-field ef- 
fect. The occurrence of new components is more pro- 
nounced in the electric field than in the magnetic part. 
Subwavelength features are present in all components but 
with very different energy levels. 
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